Griffiths phases in the strongly disordered Kondo necklace model 
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The effect of strong disorder on the one- dimensional Kondo necklace model is studied using 
a perturbative real-space renormalization group approach which becomes asymptotically exact in 
the low energy limit. The phase diagram of the model presents a random quantum critical point 
separating two phases; the random singlet phase of a quantum disordered XY chain and the random 
Kondo phase. We also consider an anisotropic version of the model and show that it maps on the 
strongly disordered transverse Ising model. The present results provide a rigorous microscopic basis 
for non-Fermi liquid behavior in disordered heavy fermions due to Griffiths phases. 



An understanding of the effects of randomness on the 
quantum critical point (QCP) of the d = 1 Kondo neck- 
lace (KN) model is relevant for the study of disor- 
dered heavy fermions systems with non-Fermi liquid be- 
havior J^, Q|. Recently a non-perturbative real space 
renormalization group (RG) was presented showing that 
weak disorder is an irrelevant perturbation near the QCP 
of a d = 1, anisotropic, pure KN model This result is 
in agreement with the generalized Harris criterion || |7j 
for irrelevance of disorder, v > 2/d, where v — 2.24 is the 
value obtained for the correlation length exponent of the 
QCP of the pure anisotropic system ||. On the other 
hand different approaches have been proposed to de- 
scribe the non-Fermi liquid behavior of disordered heavy 
fermions that rely on the relevance of disorder |J| . In 
order to settle this important point we investigate here 
the one-dimensional KN model in the case of strong dis- 
order using a generalization of a perturbative real space 
RG approach. This is the Ma-Dasgupta-Hu method Jp] 
which has been extended by Fisher j^, [l(| and others |Tl ] 
mostly for the study of random quantum spin chains and 
the random transverse- field Ising model (RTIM). In this 
Letter we obtain an exact mapping of the strongly disor- 
dered KN model into a problem of quantum spin chains. 
For the anisotropic KN model, the mapping is in the ran- 
dom transverse- field Ising model (RTIM). The general 
mapping allows to apply directly to the KN problem the 
results of extensive work done recently on random quan- 
tum spin chains §, |o[ (□], (TJ, fj||, ||, ||. This includes 
the existence of Griffiths phases with the characteristic 
singular behavior of different thermodynamic quantities 
at low temperatures. The present work provides a rigor- 
ous microscopic justification for the non-Fermi liquid be- 
havior of disordered heavy fermions due to the existence 
of a Griffiths phase in the neighborhood of a random 
QCP §. 

The one-dimensional KN model is defined by the 



Hamiltonian, 
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where and S^, fi = x,y, z are spin-1/2 Pauli matri- 
ces denoting the conduction electrons and the spins of 
the local moments, respectively. The sites i and i + 1 
are nearest-neighbors on a chain of L sites. The local 
Kondo interactions, Ji > and the hopping energies 
Wi > are uncorrelated quenched random variables with 
probability distributions, Pj(Ji) and P\y(Wi). In the 
anisotropic version of the model, (X — KN), the band 
of conduction electrons is represented just by an Ising 
term, J2±^i Wicr x icr x i + i . The full isotropic KN model 
in Eq. (|f|) will be refered from now on as the XY — KN 
to avoid confusion. 

The KN model was proposed by Doniach jjj to study 
heavy fermions and emphasizes magnetic degrees of free- 
dom neglecting charge fluctuations. It incorporates the 
essential physics of these systems which results from the 
competition between Kondo effect and magnetic order- 
ing. In the absence of disorder, the models above have 
distinct behavior. For the X — KN there is an unstable 
fixed point at a finite value of (J/W) separating an an- 
tiferromagnetic phase from a spin compensated, Kondo- 
like phase ||. For the XY — KN, any interaction J > 
gives rise to a dense Kondo state [|[ |l7| Q pl| . 

In this Letter, in order to implement the perturbative 
RG method for the KN models, we consider the conduc- 
tion electrons o~i and the local moments Si arranged in 
a chain as shown in Fig.l. Next we choose the largest 
interaction in the chain, 

n = max{Wi,Ji} (2) 

If the strongest interaction is a Kondo coupling (bond) 
between a local moment and a conduction electron, for 
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FIG. 1: The decimation processes in the case the strongest 
interaction £1 is either a) a bond or b) a hopping . 



example, il/ = Ji , the local moment 5*2 and the conduc- 
tion electron oi are decimated out yielding an effective 
hopping W between the conduction electrons <j\ and 03 
at neighboring sites which is obtained in second order 
perturbation theory (see Fig. 1). 

If the strongest interaction is a hopping term, say, 
fli = W2, the four spins 0%, Si, 03 and S3 are considered 
as a cluster. This is replaced by an effective two-spin 
cluster consisting of renormalized, local moment S and 
conduction electron a coupled by a new effective local 
Kondo interaction, J (see Fig. 1). Thus, after decimat- 
ing a strong interaction, Wi or Ji, we have an effective 
Hamiltonian with two less spin degrees of freedom and 
all couplings < Hi. 

The RG transformation gives, in the case the strongest 
interaction is a bond, an effective hopping, 



W 



Wi W 2 



and in the case it is a hopping we obtain, 
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The new Hamiltonian has exactly the same form as the 
original one, but now the system is formed by spin clus- 
ters and effective bonds. Note that the resulting flow 
equations, Eqs. (||) and (||), present a duality between 
W and J. We find that for the X ~ KN model the pa- 
rameter k = 1, so that the recursion relations for this 
model map exactly into those of the RTIM . For the 
XY - KN model, we get k = 4/V6 ~ 1.63. 

The method is implemented numerically on samples 
of sizes up to L = 2 18 and averages over 10 2 configu- 
rations. We use rectangular distributions for the local 
bonds and hopping terms. Periodic boundary conditions 
are applied. The relevant parameter is the ratio (Jo /Wo) 
of the cut-offs of the original distributions. Furthermore 
we take Wo — 1 such that Jq is taken as the control 



parameter. The dual nature of the recursion relations 
allows to locate the random QCP at Jo = 1 for any k. 

We measure the distance to the random QCP by the 
variable 



<lnJ>-<lnW^> 
var(\n J) + uar(ln W) 



(5) 



where < — > means average over quenched disorder and 
var(x) denotes the variance. Of course 5 = for Jo = 1. 

At the QCP, Jo = 1 or S = 0, the parameter k > 1 
is irrelevant and we obtain for both models a be- 
havior associated with a random singlet phase of the lo- 
cal moments. This is characterized by the fact that the 
fixed point distributions of bonds, Pj(Ji), and hoppings, 
Pw(Wi), have power-law forms, 
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with the exponents aw and aj depending on the cut-off 
il. They are given by, 



aw = otj 



1 



(7) 



The physical behavior of the thermodynamic quantities 
as a function of temperature T, which arises from such 
distributions of interactions is extensively described in 
the literature O, |2l[] . It is given by power laws with 
logarithmic corrections. 

For the anisotropic KN model, the QCP at J = 1 
separates a disordered antiferromagnetic phase (J < 1) 
from a dense Kondo compensated phase (J > 1). The 
phase for J < 1 of the XY — KN model, where Pj(Ji) 
becomes negligeable under iteration, is identified, from 
the limit Jo = 0, as a random singlet phase. In fact this 
limit corresponds to the strongly disordered XY quan- 
tum chain which, as shown by Fisher ||, exhibits a ran- 
dom singlet phase. The nature of the disordered Kondo 
compensated phase for Jo > 1 is very similar in both KN 
models and will be discussed in detail below. 

Above and below, but close to criticality, the recursion 
relations in Eqs. (j|) and (^), give rise to Griffiths phases 
with a range < \S\ < Sq. In each side of the QCP such 
phases are dominated by rare, very large clusters of the 
opposite phase. They are also characterized by power 
law behavior of the probabilities distributions, as in Eq. 
(^) , but in this case with exponents which depend on the 
distance S to the QCP. As an illustration, we show in 
Fig. 2 the exponent a,j, in the disordered phase (Jo > 1) 
as a function of the cut-off Q, as Q is reduced under 
iteration, for different distances (S < 5g) of the random 
QCP. For small 5, sufficiently close to the QCP, the 
finite intercept for (— 1/ In £7) — > is directly proportional 
to this distance <5 liTjt]. 
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FIG. 2: The exponents a j of Eq. (J6J) at the disordered dense 
Kondo phase (Jo > 1) of the X — KN model for different dis- 
tances S to the random QCP in the case of rectangular gapless 
distributions. The abscissa is the inverse of the logarithm of 
the effective cut-off. 



The range Sg of the Griffiths phases in control param- 
eter space depends on the original distributions of the in- 
teractions. For gapless distributions, i.e., those unbound 
from below, the Griffiths phases extend all over the phase 
13 1 . For bounded distributions from below, 



diagram [1* 

i.e., with a low energy gap, they have a finite extension 
around the QCP. 

The divergence of the susceptibility in the Griffiths 
phase is given by, Xo(T) oc j 11 / 2 ^ 5 )- 1 where z is the 
dynamical exponent jl^, . The behavior of the specific 
heat in this phase is given by, C V {T)/T cx r-i+V**(«) 
and for small magnetic fields H, the magnetization, 
M oc H 1 ' Zlt ^ s ' . The z exponent has the usual meaning of 
a dynamic exponent in quantum phase transitions in that 
it relates length and time scales. It is an invariant of the 
RG equations, Eqs. [|), and assumes the values z = oo 
and z = 1 at the random QCP (S = 0) and at the border 
of the Griffiths phase (\S\ = Sg), respectively p"3|, 

Fig. 3 shows the dynamic exponents z K (S) for both 
X - KN and XY - KN models in the Griffiths phase, 
at Jo > 1, for rectangular unbound distributions. Due to 
its invariance along the iteration process this exponent is 
very useful to characterize the temperature dependence 
of different thermodynamic quantities in all the Griffiths 
phase. For the original gapless distributions used to ob- 
tain z K (5) in this figure, this phase extends for all 6 > 
jl3], Oj. Note from Fig. 3 that for sufficiently small S, 
* K (<5pcx<S|,p. 

Starting the iteration with a distribution Pj(Ji) with 
a gap, in the disordered Kondo phase (Jo > 1) we obtain 
for both the X-KN and XY- KN, besides the Griffiths 
phase for S < Sg, a strongly disordered random Kondo 



FIG. 3: The inverse of the dynamic exponents z for the 
X — KN and XY — KN models as function of the distances 
to the random QCP (S > 0). They are obtained iterating 
numerically the recursion relations starting from gapless rect- 
angular distributions. The dynamic exponents determine the 
singularities of the thermodynamic quantities with decreasing 
temperature in the Griffiths phases (see text). 



phase (RKP) for 8 > 8q- This phase consists essentially 
of a collection of isolated Kondo singlets with a distribu- 
tion of excitation energies. It is natural then to describe 
such phase using a distribution of Kondo temperatures 
§■ 

In Fig. 4 we show the uniform susceptibility as a func- 
tion of temperature, x(T), for the X — KN model, cal- 
culated far away from the random QCP, deep in the dis- 
ordered Kondo phase for two initial distributions Pj{Ji). 
When the original distribution of Kondo couplings Jj has 
a finite gap Aj at low energies, the susceptibility (solid 
line) does not diverge at low T as shown in the figure. In 
this case the weakly disordered Griffiths phase extends to 
5 < Sg — 0.358 which corresponds to a gap in the origi- 
nal Pj(Ji) of Aj < 0.156. The susceptibility shown is for 
S = 0.44 > Sg, in the RKP. It is clear that in this case 
the fixed point distribution of gaps or Kondo tempera- 
tures does not become sufficiently singular in the RKP 
to yield a diverging susceptibility. In the same figure 
we show the temperature dependent susceptibility for an 
initial unbound distribution of interactions, i.e., a rect- 
angular distribution with no gap. The susceptibility now 
diverges as T — * 0. Since in this case of gapless distribu- 
tions, the Griffiths phase extends all over the disordered 
Kondo region, S > 0, the divergence of x(^) is clearly 
associated with this Griffiths phase. Similar behavior is 
found for the XY — KN model. 

In summary, we have studied the effects of strong dis- 
order in the KN model for heavy fermion systems using 
a perturbative RG approach. Differently from the case 
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The dense Kondo phase for Jq > I (5 > 0) is similar 
in both KN models studied. For the XY - KN with 
Jo < 1 (5 < 0) there is a random singlet phase coexist- 
ing for some range of <5, which depends on the original 
distributions, with a weakly disordered Griffiths phase 
pp| . The exact mapping of the recursion relations for 
the X — KN on those of the RTIM allows to carry on 
the extensive results obtained on the latter to the present 
former model. It enables us to show unambiguously the 
existence of non-Fermi liquid behavior in strongly disor- 
dered heavy fermions associated with Griffiths phases. 

The authors are grateful to the Brazilian agencies 
FAPERJ, CNPq and CAPES for financial support. 



FIG. 4: The uniform susceptibility of the X — KN model as a 
function of temperature for two initial distributions of bonds 
and the same distance S — 0.44 of the random QCP. For 
a gapped distribution (solid line), for which S > 5g ~ 0.36, 
and the system is in the RKP. For a gapless distribution 
(dashed line) such that the system is in the Griffiths phase. 
Temperature is in units of Wo- 



of weak disorder || , strong disorder is a relevant pertur- 
bation and gives rise to Griffiths phases which for distri- 
butions with a gap extend over finite regions of the phase 
diagram, above and below the QCP. For Jo > 1, in the 
strongly disordered regime above the Griffiths phase, we 
find a random Kondo phase which consists of a collec- 
tion of isolated singlets with random excitation energies. 
It is natural to describe such a phase by a model with 
a distribution of Kondo temperatures The suscep- 
tibility however does not diverge at low temperatures 
in the RKP. We observe a diverging susceptibility in 
the case of lower unbound distributions where the Grif- 
fiths phase extends all over the disordered region. Conse- 
quently our results indicate that a diverging x(T) is due 
to a Griffiths phase. It can not be obtained in a model 
of isolated Kondo singlets with random Kondo tempera- 
tures. A fixed point distribution Pj(Ji) of the power law 
type, that gives rise to a diverging x(P)i is incompatible 
with the existence of isolated singlets. In the Griffiths 
phases, the temperature dependence of the thermody- 
namic quantities is expressed in terms of the dynamic 
exponent z. We have obtained this exponent for both 
KN models as a function of the distance to the QCP for 
gapless distributions. 

At the random QCP of both the X - KN and the 
XY — KN there are random singlet phases character- 
ized by thermodynamic quantities with power law be- 
havior, arising from distributions of interactions like in 
Eqs. (|^) and (0), with weakly temperature dependent 
exponents. For the disordered X — KN model there is 
a disordered antiferromagnetic phase for Jq < 1 (6 < 0). 
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